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We study transport of local magnetization in a Heisenberg spin-1/2 chain at zero temperature. 
The system is initially prepared in a highly excited pure state far from equilibrium and its evolution is 
analyzed via exact diagonalization. Integrable and non-integrable regimes are obtained by adjusting 
the parameters of the Hamiltonian, which allows for the comparison of transport behaviors in both 
limits. In the presence of nearest neighbor interactions only, the transport behavior in the integrable 
clean system contrasts with the chaotic chain with on-site defects, oscillations in the first suggesting 
ballistic transport and a fast decay in the latter indicating diffusive transport. The results for 
O ■ a non-integrable system with frustration are less conclusive, similarities with the integrable chain 

being verified. We also show how methods of quantum control may be applied to chaotic systems 
I ) \ to induce a desired transport behavior, such as that of an integrable system. 
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(p ; I. INTRODUCTION 

^ The understanding of the behavior of quantum many-body systems far from equilibrium is still far from satisfactory. 

• Under debate, for example, are the conditions that lead to the thermalization of isolated quantum systems [11, 0, 0| 
^ or to ballistic and diffusive transport. The distinction between ballistic and diffusive transport in linear response 
(_i ■ theory usually relies on the computation of the Drude weight, which is finite in the first case 0,0], but alternatives 
approaches have also been considered. In studies of thermal transport, for instance, the derivation of the Fourier law 
H from a microscopic foundation has been pursued by applying the Hilbert space average method @ or by numerically 
' simulating open systems coupled to heat reservoirs [7|, |8| . 

In addition to fundamental questions, the subject of quantum transport received increased attention after experi- 
mental observations of unusual excess conductivity for heat d, [l^, [ll|, and magnetization [l^ in low-dimensional 
Q ' magnetic compounds. Advances in optical lattices are also very promising, since these highly controllable systems 
allow for the investigation of the dynamics of strongly correlated quantum systems for long times. Ultracold atom 
gases have been realized in optical lattices [l^ and magnetic systems may soon be implemented [H, [l3j ■ Magnetic 
^ , compounds are well described by Heisenberg spin-1/2 systems, which in part explains the proliferation of articles 
(^s^ ' analyzing transport in these systems. Under certain conditions, this model may also be mapped into spinless fermion 
. systems via a Jordan- Wigner transformation fisj or onto systems of hard-core bosons 0| ■ 

' A special feature of the one-dimensional Heisenberg spin-1/2 model in the absence of defects and with nearest- 
neighbor interactions only is that the model is integrable and solved with the Bethe Ansatz method 19]. The 
addition of impurities [20, [2l| or next-nearest neighbor interactions [l^, [2^ to the chain may lead to the onset of 
^ quantum chaos. By varying the parameters of the Hamiltonian we may then analyze transport behavior in both 
0^ . regimes. At finite temperature, in the gapless phase, o ne g enerally associates integrability with ballistic transport 
' and chaoticity with diffusion [^, [23, [H, [la [13, [H, [3, [30| . In the gapped phase, diffusive behavior is found also 
J> in integrable systems [Sll. Is^] . However, there are still open questions, especially about the transport in integrable 
systems at finite [sol |32|] and infinite temperatures [ssj , and in nonintegrable systems at low temperatures [34t ISSl . [36j 
. , - such as the possibility of ballistic transport in frustrated chains [13, [lo, [H, [13, [H, [H, [10, [4l| . At zero temperature, 
' within linear response theory the correspondence 'gapless phase - ballistic transport' and 'gapped phase - diffusive 
■ ■ '■ behavior' is well established [sil. [4^. but a new discussion has emerged concerning the problem of quench dynamics 
in the case of initial states very far from equilibrium. The latter is the focus of the present work. 

Here, we compare transport of local magnetization in an isolated Heisenberg spin-1/2 chain with free boundaries 
in both integrable and non-integrable regimes. A highly excited initial state corresponding to spins pointing up in 
the first half of the chain and pointing down in the other half is considered. We assume, as in Ref. [29], that an 
exponential decay of the magnetization to equilibrium indicates diffusive transport and that a bouncing behavior 
indicates ballistic transport. Two chaotic chains are considered: clean with frustration in the three directions and 
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disordered with a single defect in the middle of the chain. Exact diagonalization is used, which limits the analysis to 
small system sizes. An exponential relaxation to equilibrium is verified for the disordered system, while the frustrated 
chain shows oscillations similar to those seen in the integrable system [i^ . It has been suggested that such oscillations 
may not be related to ballistic transport [3^ , but are possibly related to the quantum phase transition of the isotropic 
chain [13, El]- We show that even though reduced, the oscillations persist also for anisotropic frustrated systems in 
the gapless phase. We argue that the similarities between frustrated and integrable chain result from the mixing of 
different symmetry sectors, which does not occur in the dynamics of the disordered system. In the gapped phase, 
localization of local magnetization is observed for all systems, as expected. 

In the second part of this work, we discuss how methods of quantum control may be used to induce a transition from 
diffusive to ballistic transport. This is accomplished by applying a sequence of control operations which eliminate 
unwanted terms of the Hamiltonian [4^ 113, . The system studied here is subjected to a static magnetic field in 
the z direction; all sites having the same Zeeman splitting, except the defect. By applying a sequence of very strong 
magnetic fields that frequently rotate all spins perpendicularly to the z direction, we show that the effects of the 
on-site disorder may be averaged out and the transport behavior of the integrable chain recovered. A sequence that 
eliminates the effects of next-nearest neighbor interactions has also been proposed in ^is^ . 

The paper is organized as follows. Sec. II describes the model and the quantity used to characterize chaos. Sec. Ill 
shows the results for local magnetization in the gapless and gapped phase for both integrable and non-integrable 
regimes. Sec. IV compares the transport behavior of the disordered chain in the absence and presence of a control 
sequence. Discussions and concluding remarks are presented in Sec. V. 



II. SYSTEM MODEL AND CHARACTERIZATION 



We consider a one-dimensional Heisenberg spin- 1/2 system with open boundary conditions described by the Hamil- 
tonian, 



H — Hz + Hnn + Hnnn 

L L-l L-2 

n—1 n—1 n—1 

where h — 1, L \s. the number of sites, and S'i^'^'^^' are the spin operators at site n in the three directions. The 
parameter Un is the Zeeman splitting of spin n as determined by a static magnetic field in the z direction. The system 
is clean when all sites have the same energy splitting w„ = w, and it is disordered when defects characterized by 
different energy splittings w„ = w -I- d„ are present. J and J' are the exchange coupling strengths of nearest-neighbor 
(NN) and next-nearest-neighbor (NNN) couplings, respectively, and are assumed to be constant; A is the anisotropy 
associated with the Ising interaction S^^S^^^^. At zero temperature, — 1 < A < 1 corresponds to the gapless phase 
and |A| > 1 to the gapped phase. We set J, J', A > 0. In the following, a spin up may sometimes be referred to as 
an excitation. 

The total spin operator in the z direction, — X]n=i ^rn conserved, so the matrix H is composed of independent 
blocks, each of dimension N — (j^) — i!/[(L — Af)!M!], where M is the total number of spins pointing up. We consider 
L even and study the largest subspace, = 0, which corresponds to M = L/2 (from the perspective of interacting 
spinless fermions, we are at half-filling). The calculations are performed in the basis consisting of eigenvectors of . 

In quantum systems, integrable and non-integrable regimes may be identified by analyzing the distribution of 
spacings sp between neighboring energy levels [49i [soj . Quantum levels of integrable systems tend to cluster and are 
not prohibited from crossing, the typical distribution is Poissonian: Pp{sp) = exp(— sp). In contrast, non-integrable 
systems show levels that are correlated and crossings are avoided, the level statistics is given by the Wigner-Dyson 
distribution. The form of the Wigner-Dyson distribution depends on the symmetry properties of the Hamiltonian. 
In the case of systems with time reversal invariance we have: Pwd{sp) = (7rsp/2) exp(— 7rsp^/4). The Heisenberg 
model with a magnetic field, as in ([T]), does not commute with the conventional time-reversal operator, but the 
Hamiltonian elements are real and the distribution is still given by Pwd{sp) pollsH. 

A meaningful level spacing distribution requires the separation of different symmetry sectors. If energies from 
different subspaces are mixed, we may obtain a Poisson distribution even if the system is chaotic. The main feature 
of chaoticity, level repulsion, is therefore lost. In addition to 5^, H ([Ij conserves also parity, when the system is 
clean, and total spin 5^ ~ (X]n=i ^nY, when the system is clean and isotropic, A = 1. A particular feature of the 



= Q subspace is its invariance also under a 7r-rotation of all spins 51| . This symmetry is taken into account when 
S"^ is considered, but it remains for an anisotropic system [23,]. In the case of a disordered chain with NN interactions 
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FIG. 1: (Color online.) Onset of chaos in a Heisenberg chain with L — 12. Left panel: Disordered chain with NN interactions 
only, as described hy H ^ with J' — and ~ 0.5 J. Both A = 1 and A — 0.5 lead to distributions similar to Pwd- 

Dotted region - Pp; Horizontal lines - Pwd- Middle panel: Frustrated chain, as described hy H ^ with J' = J. The chaotic 
system with A = 1 and d„ = gives a Poisson distribution when and reflection symmetries are mixed. A Wigner-Dyson 
distribution is recovered when both symmetries are broken by considering A — 0.8 and (11 = 0.5J. Right panel: Histogram 
of eigenstates of H ^ with J' = J and A = 1 according to their 5*^ subspaces. Empty bars: all eigenstates; filled bars: 
eigenstates j that contribute to with probability amplitude |aj,in| > 10~^. 



only, a single defect placed in the middle of the chain, dL/2+1 7^ 0, may lead to the onset of chaos |21|, as shown 
in the left panel of Fig. [TJ A clean system with J = J' and A = 1 is also chaotic [l^l, but it conserves and 
parity. If these symmetries are not taken into account, a misleading Poisson distribution is obtained, as shown in the 
middle panel of Fig. [TJ How these states distribute among the various S'^ subspaces is given in the right panel for 
L = 12; moreover both parities are present, from the N — 924 eigenstates, 472 have odd parity. The symmetry 
may be broken by choosing A ^ 0, but here spin reversal and translation invariance remain. The three symmetries 
are guaranteed to be broken by placing a defect in one of the edg es of the chain Jl,, ,5^] , an alternative that does not 
affect the integrability of the clean system with NN-coupling [531 ■ By combining both, defect and anisotropy, a clear 
Wigner-Dyson distribution is reached for the system with NNN interactions, as shown in the middle panel. 



III. TRANSPORT BEHAVIOR: LOCAL MAGNETIZATION 



We study transport of local magnetization as defined by 



L/2 



M(t)^(vl/(t)|^^;j*(t)), 



(2) 



where is the state of the system at instant t. 

The initial state |5'(0)) considered here is the basis vector {(pin) of corresponding to spins pointing up in the 
first half of the chain and pointing down in the other half. The eigenvectors jV'} of 77 ([T]) in the basis {(f)) are given by 
= aj^k\4'k) with eigenvalues Ej, so the local magnetization may be written as 



Mit) 



1 ^ 



L/2 \ / ^ ^ 
in=l / \fc=l p=l g=p+l 



cos[{Ep - Eq)t] 



(3) 



where = —1 (+1) if spin n of state is pointing down (up). 

The scenario is that of a quench dynamics, where the system is initially prepared in the eigenstate |5'(0)) = | ff 
. . . II ... II) of Hz and then suddenly interactions are turned on and the system dynamics becomes dictated by 

Integrable vs. disordered chain. The clean system with NN interactions only, as described by H ^ with d„ — 
and J' = 0, is an integrable model. The dynamics of local magnetization for different values of A is shown in the top 
panels of Fig. [2| When A > 1, few states are in resonance with \(f>in) and therefore localization occurs (right panel). 
In the gapless phase, the bouncing behavior of M{t) suggests ballistic transport [2^. However, it calls attention the 
reduction of the oscillations amplitudes as A decreases (left panel), which is probably associated with the breaking 
of the symmetry. 
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FIG. 2: (Color online.) Transport of local magnetization in a Heisenberg chain with NN interactions described by -ff HI with 
J' — and L — 12. The value of u is irrelevant for the dynamics and A is indicated. Top panels: M{t) for the clean integrable 
system with d„ = 0. Bottom panels: M{t) for the chaotic disordered system with dL/2+i = 0.5 J. 



The addition of on-site disorder to the chain with NN interactions breaks S*^ symmetry and parity. A single defect 
d'L/2+i = 0.5 J leads to the onset of chaos [S^l and to a fast decay of the local magnetization to equilibrium for 
0.4 < A < 0.8, as seen in the bottom left panel of Fig[2l As A increases - bottom right panel - defect and anisotropy 
reduce resonances with the initial state and localization eventually takes place. 

Integrable vs. frustrated chain. H ([1]) with rf„ = and J = J' describes a chaotic system. Instead of diffusive 
transport, however, the plots in Fig. [3] indicate the following, (i) Localization occurs again, as expected, when A > 1. 
(ii) When A ~ 1, the evolution of local magnetization shows a bouncing behavior that is similar to that observed for 
the integrable system. The oscillations for J — J' are faster than for J' = 0, but their amplitudes are of the same 
order (cf. top right panels of Figs. [2] and [3]). The bottom left panel suggests that these oscillations are not artifacts 
of the small size of the chain, since their amplitudes do not appear to decrease as the system size increases, although 
simulations for larger systems are necessary to settle this point [5^ . (iii) According to the top left panel of Fig. [3l the 
reduction of oscillations for the system with NNN interactions associated with the decrease of anisotropy, A < 0.8, 
is more significant than in the integrable case (cf. top left panels of Figs. [Hand [3]). Notice that the curve for the 
frustrated chain with A = 0.2 is hardly noticeable due to the absence of oscillations and the fast decay of M{t) to 
equilibrium. 

We interpret observations (ii) and (iii) as follows. State |^'(t)) is a linear combination of eigenstates of H belonging 
to a single S*^ sector, but no restriction exists with respect to S"^, global tt rotation, and parity. Similarly to what 
is observed for the level spacing distribution (middle panel of Fig. [T]), but now from a dynamical point of view, the 
effects of level repulsion may be missed by the local magnetization if those symmetries are present. The mixing of 
symmetry sectors implies that a large fraction of (almost) degenerate states contribute to M{t), such level crossings 
must reduce the phase randomization that leads to a fast decay of the local magnetization ([3]). For instance, when 
i — 12 and J' — J, there are 429 eigenstates \ipj) participating in |5'(i)} with probability amplitude |aj,i„| > 10~^. 
From 400 states (the border ones being neglected), there are 29 unfolded spacings between neighboring levels with 
sp < 0.1. These 429 eigenstates are distributed among the S'^ subspaces with quantum number S between and 
21 as shown in the right panel of Fig. [U and both parities are found, 217 states having odd parity. This is to be 
contrasted with the disordered system, where is the only conserved quantity. M{t) captures well the transition to 
chaos and diffusive transport becomes evident. For L = 12, J' = 0,A = 1, and dj = 0.5 J, there are 751 eigenstates 
IV'j) participating in \'^{t)) with probability amplitude |aj_i„| > 10^*^. From the 720 states in the bulk, only 5 lead to 
sp < 0.1, a percentage 10 times smaller than in the frustrated chain. 

For the frustrated chain, as A decreases from 1 and the S*^ symmetry is broken, transport of local magnetization 
becomes closer to a diffusive behavior. However, since parity and invariance under a global tt rotation are still present, 
the decay of M is not as abrupt as in the disordered system. Better agreement is verified if these symmetries are also 
broken by adding, for example, a defect to one of the borders of the chain, as in the bottom right panel of Fig. [S] 
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FIG. 3: (Color online.) Transport of local magnetization in a Heisenberg chain with NNN interactions described hy H ^ with 
J' — J and L = 12. Top panels: M{t) for the frustrated system with different values of A. Left bottom panel: Comparison of 
M{t) for different system sizes, A = 1. Right bottom panel: Approach to diffusive transport as symmetries are broken. Present 
symmetries: and parity when A = 1 ; parity when A = 0.8; none of the two when A = 0.8 and = 0.5J. 



This behavior should again be compared with the middle panel of Fig. (TJ where P{sp) approaches a Wigner-Dyson 
distribution when the symmetries are broken. 



IV. MANIPULATION OF TRANSPORT BEHAVIOR VIA QUANTUM CONTROL 

We have recently discussed the use of quantum control methods, such as dynamical decoupling (DD), to induce 
a desired transport behavior [i^. DD techniques consist of sequences of external control operations. The goal is 
to average out the effects of unwanted terms and reshape the Hamiltonian in order to achieve a desired dynamics. 
These methods were first applied in Nuclear Magnetic Resonance (NMR) spectroscopy, where the control operations 
correspond to very strong magnetic fields (pulses) able to rotate the spins and time-reverse the system evolution (45l . 
l46j . More recently their applications have been extended to more general systems, especially in the context of quantum 
information, where the purpose is not only control of internal interactions [56j , but also of external couplings [47i [s^ . 

Pulse sequences have been developed and used for several different systems. Sequences to eliminate NNN interactions 
or on-site disorder from H ((!]) were proposed in Ref. [i^. Consider, for example, the chaotic system with on-site 
disorder, be it a single defect or various impurities spread through the chain In order to eliminate the effects of 
the disorder and recover the transport behavior of the integrable system, we apply sequences of instantaneous pulses 
given by 



= exp ^ S"^ = exp(-i7rS'^), 



\ n=l / 

which rotate all spins 180° around the x direction. The pulses at instants tj-i and tj are separated by a time interval 
T — tj — tj^i of free evolution. They alternate the sign of the one-body terms, but do not change the two-body 
interaction terms. The propagator at time t2p — pTc, where p £ N and Tc = 2t is the cycle time, approaches that of 
an integrable clean system as t — > 

U{pTc) eyi^[-iHNNpTc]. 
This result is achieved as follows. Consider the notation 



U+ = exp[-z(i/2 + Hnn)t] = exp[-ii7ir], 
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f/_ = P^U{tj,tj_i)P^ = P^{P,Pl)U{tj,tj_i)P, 



-1 exp 



(« 



-I e 



(H, + HNN)e 



where 1 is the identity operator. The propagator at t2p can then be written as 



UipT,) = P,U{hp,t2p-i)P,...P,U{t2,h)P,U{h,0) 
= U^U+...U^U+^ exp[-ii?pTj. 



(4) 



The average Hamiltonian, H 
The first three dominant terms are 



Sfc^o^''^'' obtained by using the Baker-Campbell-HausdorfF expansion [isl. |46 



C 

As r ^ 0, the average Hamiltonian then approaches Hj^n and the transport of local magnetization coincides with 
the one obtained for the integrable chain. For small values of r, the agreement may hold for relatively long times, 
as shown in Fig. 3) In real scenarios, however, we need to take into account several pulse imperfections and the 
limitation on pulse separations, which leads to the accumulation of residual averaging errors. Different approaches 
exist to handle these problems (see Ref. [HI and references therein). 




400 



FIG. 4: (Color online.) Transport of local magnetization in a Heisenberg chain described hy H ^ with J' = 0, A = 1, and 
L = 12. The curve showing a fast decay of M{t) corresponds to a disordered system with — Q.35J in the absence of pulses. 
Bouncing curves represent the clean chain (dashed line) and the disordered system (full line) subjected to the DD sequence 
Q. Data is acquired after every Tc = 2t, where t = J~^. 



CONCLUSION 



We compared the transport behavior of local magnetization in integrable and chaotic one-dimensional many-body 
systems in the case of an initial state far from equilibrium and studied quantum control methods as possible tools 
to manipulate transport behavior. An isolated Heisenberg spin-1/2 chain was the study-case. Chaos was induced by 
adding on-site disorder or frustration to the system. 

We found that in the gapped phase, local magnetization is localized for all systems. In the gapless phase, diffusive 
transport is evident in the disordered chain, but the frustrated chain shows a behavior similar to the integrable model. 
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especially when closely isotropic. We argued that this is a consequence of the mixing of different symmetry sectors, 
which occurs in the dynamics of both the frustrated and the integrable system. The combination of different subspaces 
obfuscates the effects of level crossing, which is a main feature of chaotic systems, and slows down the decay of local 
magnetization to equilibrium. These results indicate that a one to one correspondence between quantum chaos and 
diffusive transport docs not necessarily hold. To be able to view transport as a dynamical signature of chaos, we need 
to consider quantities that reside in a single symmetry sector. Local magnetization, for instance, is such a quantity 
when parity and total spin are not conserved, as in the disordered system, since it resides in a single subspace of total 
spin in the z direction. 

We also discussed a sequence of control operations that may be applied to a disordered system in order to remove the 
effects of the defects and recover ballistic transport. The possibility to manipulate transport behavior via quantum 
control methods may find practical applications, such as the design of devices with controllable conductivity, the 
mitigation of the effects of local heating, or the transport of information in quantum computers [sol . [60j . Real 
systems where these ideas might be tested include crystals of fluorapatite [s^ [g^l, as studied in NMR, or optical 
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